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Abstract 



We analyze the stability properties of the purely magnetic, static solutions to the Einstein- 
Yang-Mills equations with cosmological constant. It is shown that all three classes of so- 
lutions found in a recent study are unstable under spherical perturbations. Specifically, we 
argue that the configurations have n unstable modes in each parity sector, where n is the 
number of nodes of the magnetic Yang-Mills amplitude of the background solution. 

The "sphaleron-like" instabilities (odd parity modes) decouple from the gravitational 
perturbations. They are obtained from a regular Schrodinger equation after a supersymmet- 
ric transformation. 

The body of the work is devoted to the fluctuations with even parity. The main diffi- 
culty arises because the Schwarzschild gauge - which is usually imposed to eliminate the 
gravitational perturbations from the Yang-Mills equation - is not regular for solutions with 
compact spatial topology. In order to overcome this problem, we derive a gauge invariant 
formalism by virtue of which the unphysical (gauge) modes can be isolated. 



1 On leave of absence from Tbilisi Mathematical Institute, 380093 Tbilisi, Georgia 



1 Introduction 



A previous numerical analysis |T] has revealed the following features of the purely magnetic 
solutions to the spherically symmetric Einstein- Yang-Mills (EYM) equations with positive 
cosmological constant A: The static configurations fall into three classes, where the solutions 
in each class are characterized by the values of A and the number of zeroes, n, of the magnetic 
YM amplitude w. For sufficiently small values, A < A crit (n), the solutions asymptotically 
approach the deSitter geometry and can be viewed as Bartnik-McKinnon (BK) solitons || 
surrounded by a cosmological horizon. For a set of values A reg (n) exceeding A crit (n), the 
configurations have the topology MAS 3 , where the ground state (n = 1) is the Einstein 
Universe with constant YM energy density. Between A cr i t (n) and A reg (n) there exists a 
discrete family of "bag of gold" solutions, that is, of configurations with horizon and "finite 
size" . 

This paper is devoted to the stability properties of the solutions found in ]1| with respect 
to spherically symmetric fluctuations. The metric g of a spherically symmetric spacetime is 
described in terms of the metric g of a 2-dimensional pseudo-Riemannian manifold M and 
a function R on M, such that g = g + R 2 dQ 2 . 

The stability analysis of solutions for which R has no critical points is considerably 
simplified by the fact that R can be used as a coordinate on M. Moreover, it turns out to be 
sufficient to analyze the perturbations in the Schwarzschild gauge, 5R = 0. The linearized 
Einstein equations then assume the form of constraint equations, by virtue of which one can 
express the metric perturbations 5g in terms of the fluctuations of the matter fields. In this 
way, one always ends up with a set of fluctuation equations for the matter fields in standard 
form 0. For the BK solitons and the corresponding black hole solutions 0, [|], ||, the 
set of pulsation equations therefore reduces to a single regular Schrodinger equation for the 
YM amplitude w. The number of bound states of this equation eventually determines the 
number of unstable modes M, || (in the even parity sector). 

In the present paper we are mainly interested in the stability properties of solutions for 
which R does have a local extremum. In this case, the Schwarzschild gauge, 5R = 0, is 
not everywhere regular and the above procedure can therefore not be readily adopted. In 
fact, for 5R ^ 0, not all of the linearized Einstein equations have the form of constraints. 
Hence, one obtains a coupled system of pulsation equations, including both matter and 
metric perturbations. Although this system is self-adjoint, the differential operator is not 
hyperbolic, reflecting the fact that the equations contain unphysical degrees of freedom (gauge 
modes). 

Studying the gauge dependence reveals that the fluctuation equations actually contain 
only one degree of freedom which is of physical nature. In fact, it turns out to be possi- 
ble to derive a pulsation equation for the gauge invariant quantity ( associated with the 
physical degree of freedom. The procedure generalizes the method described above for the 
Schwarzschild gauge in a gauge invariant way. 

Although the pulsation equation for ( is a standard Schrodinger equation, it has the flaw 
that the potential is unbounded. In fact, it turns out that the equator (i.e., the critical point 
of R) is a regular singular point of the Sturm-Liouville equation for (. A closer look reveals 
that the so-called limit point case occurs. "Weyl's alternative" (see, e.g. ||) then implies 
that there exists exactly one essentially self-adjoint realization of the pulsation operator, 
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and that £ must vanish at the equator for this realization. However, there exist analytic 
solutions which do not vanish at the equator and, nevertheless, give rise to perfectly regular 
metric fluctuations. 

In order to reconstruct the metric and matter perturbations from the gauge invariant 
quantity (, we consider a class of gauges for which all fluctuations can be obtained by solving 
a set of ordinary differential equations. A particular gauge of this type is the Schwarzschild 
gauge, which is, however, not regular at the equator. A globally regular gauge is obtained 
by considering conformal perturbations of the spatial part of the metric. In this way one 
finally obtains a regular set of metric and matter fluctuations from each regular solution of 
the gauge invariant pulsation equation. 

The entire procedure is nicely illustrated for the lowest (n = 1) compact solution, since 
all steps can be performed analytically in this case. The gauge invariant equation has one 
bound state, corresponding to an unstable mode. A numerical analysis reveals that the higher 
compact solutions (n > 1) have, as expected, exactly n unstable (gravitational) modes. 

The stability analysis described above applies to the family of compact solutions, para- 
metrized by n and A reg (n). The symmetry properties of these configurations with respect to 
reflections at the equator simplifies the discussion of the regular singular point. A careful 
numerical analysis shows that the solutions with horizon (bag of gold solutions) share the 
stability properties of the compact solutions. 

The stability properties of purely magnetic, spherically symmetric, static solutions to the 
EYM equations are determined by two, completely independent sets of perturbations. The 
investigation discussed above applies to the even-parity fluctuations, which are also called 
"gravitational" perturbations, because they have no flat spacetime analogues. Considering 
fluctuations of the non-magnetic parts of the YM potential yields a second, orthogonal set of 
perturbations with odd parity. For the BK solitons and black holes, these "sphaleron-like" 
perturbations give rise to additional n unstable modes (see JTDfl, |TT[ and references therein). 



Adopting the method developed in JO] , we show in the last part of this paper that this result 
remains true for both the compact solutions and the solutions with horizon found in [T|]. 

This article is organized as follows: In the second section we give the basic equations for 
spherically symmetric EYM models and recall the most important features of the compact 
solutions to these equations. 

The third section is devoted to the derivation of the pulsation equations. We start by 
linearizing the YM equations within the purely magnetic ansatz. We then argue that - 
by virtue of the Bianchi identity - the (00) and (01) components of the Einstein equations 
give only rise to one independent linearized equation, which has the form of a constraint. 
Finally, we linearize the partial traces (with respect to g and the metric of S 2 ) of the 
Einstein equations, which yields two pulsation equations for the metric fluctuations. The 
four equations obtained in this way can be combined into a self-adjoint system. However, 
since no particular gauge was fixed, the system still contains gauge modes, as is also indicated 
by the form of the differential operator. 

In the forth section we argue that there exists only one gauge invariant quantity, (, say, in 
the even parity sector. Unfortunately, the Sturm-Liouville equation for ( has - in addition to 
the origin - a regular singular point at the equator. Since both branches of the fundamental 
system are analytic in the vicinity of the equator, it is not sufficient to consider only the 
solutions belonging to the unique self-adjoint extension of the fluctuation operator (i.e., the 
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solutions for which ( vanishes at the equator). 

The reconstruction of the perturbations from the gauge invariant quantity is presented 
in the fifth section. Using the globally regular conformal gauge and the general pulsation 
equations derived in the third section, we give explicit, regular integral expressions for the 
metric and the matter fluctuations in terms of (. We also explain how to use the residual 
gauge freedom in order to fix the integration functions. The results of the numerical analysis 
are presented at the end of this section. 

An analytic discussion of the stability properties of the compact ground state solution 
in = 1, A = A reg (l)) is presented in the sixth section. It turns out that the unstable mode 
corresponds to spatially constant, conformal perturbations of the spacetime metric. This 
"breathing" mode is also obtained by an alternative approach, which takes advantage of the 



special geometry of the ground state, in order to introduce the Bardeen potentials |T2[ used 
in cosmological perturbation theory (see, e.g., |0§). We conclude this section by deriving 
a "dual" Schrodinger equation for the perturbations of the ground state. In contrast to 
the equation for (, the new pulsation equation for the "supersymmetric partner" of ( is 
completely regular at the equator. 

In the last section of this paper we consider the perturbations with odd parity. Again, 
the fluctuation equations can be reduced to a single Schrodinger type equation. Since the 
potential in this equation is unbounded at the zeroes of the YM amplitude w, we use the 
residual gauge freedom to perform a supersymmetric transformation. In this way, we obtain 
a pulsation equation with everywhere regular potential. We then demonstrate that the 
zero energy solution to this equation has the same number of nodes as w. This eventually 
proves that both the compact and the bag of gold solutions have n unstable modes in the 
"sphaleron" sector. 

2 Preliminaries 

In this paper we consider the question of linear stability for the solutions of the EYM 
equations with cosmological constant presented in |J. Before we do so, we shall briefly 
recall the basic notions and equations used in |IJ . Since we are particularly interested in the 
globally regular, compact configurations, we also recall their basic features. 

2.1 Static Equations 

The metric of a spherically symmetric manifold (M, g) can be written in the form 

g = g + R 2 g, (1) 

where g is the metric on the pseudo-Riemannian manifold M and g denotes the standard 
metric on S 2 . The Einstein tensor for the metric ([I]) becomes (see, e.g. ||) 

(2) 
(3) 
(4) 
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where the quantities with a tilde refer to (M, g) and those with a hat to (S 2 ,g). (We use small 
and capital Latin letters for indices on (M,g) and (S 2 ,g), respectively; a, b, c = 0, 1 and 
A, B, C = 2, 3.) Without loss of generality, we shall often use the diagonal parametrization 



g = -e 2a ^dt 2 + e 2b ^dp 2 



(5) 



of the metric on M. With respect to this, the d'Alembertian of a function (e.g., R) and the 
Ricci scalar on (M,g) become (with R = dR/dt, R' = dR/dp) 
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(7) 



respectively. 

A purely magnetic, spherically symmetric 577(2) YM gauge potential is parametrized in 
terms of a scalar function w : M — > M (see, e.g. 0) 



A = (w — 1) [ftp d-d — % sm^d dtp] 
The stress-energy tensor for A has the components 
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where g 2 is the gauge coupling constant and V(w) = (1 — w 2 ) 2 . 
Using the notation 

R = exp(ji) , 

and the parametrization (|I) of the metric g, the static field equations are 



(8) 

(9) 
(10) 
(11) 

(12) 

(13) 
(14) 
(15) 

AR 2 — (16) 

Note that eqs. (U), flTJj) and (|TJ) are the |(00 + 11), |(00 - 11) and §(00 + 11 - 22 - 
33) components of the Einstein equations. We also recall that the dimension-full coupling 
constant k = 8nG/g 2 can be absorbed by using the dimensionless quantities R/yfn, p/y^ 
and Ak (see [|T[). Hence, we use k = 1 throughout. (This corresponds to Johnstone Stoney 
units, already used in 1881 fL"4|.) 
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2.2 Compact Solutions 

In the first six sections of this paper we are mainly interested in the stability of the com- 
pact, regular solutions. These are characterized in terms of the node number n of the YM 
amplitude w and the value of the cosmological constant, A = A rcg (n). The lowest solution 
belongs to A = 3/ (2k) and can be given in closed form, 

a = , 6 = 0, Rj \l~k~ = sin(p/ ^/k) , 

w = cos(p/v/k), with A = 3/(2«). (17) 

The solution describes the static Einstein Universe, M x S 3 , with T 00 = \g~ 2 -^ and T u = 
T22 = T 33 = \g~ 2 -^- (For later use we also note that - for k — 1 - we have p! = cot p, 
V(w) = sin 4 p, and V, w / (4R 2 ) = \3w = — cos p.) 

We recall that there are no solutions of the static equations for which // has more than 
one zero. (The first non- vanishing derivative of \i is negative for every p with p!(p) = 0; see 
the expansions below.) The unique value p e with p\p e ) = (and R(p e ) 7^ 0) is therefore 
called the equator. There are two families of solutions which are analytical in the vicinity 
of the equator. According to the behavior of the YM amplitude w, these solutions will be 
called odd and even. 

For the odd configurations, w (p e ) = 0, w'(p e ) 7^ 0, one finds with x = p — p e (and k = 1) 
fi' = -(^) 2 x + 0(x 3 ), w' = w' e + 0(x 2 ), 

Re 

a = a e + 0(x 2 ) , b = b e + 0{x 2 ) , (18) 
where the parameters w' e and R e are subject to the constraint equation 

KRt - R 2 + \ = e' 2b (R e w' e ) 2 . 

The even solutions, for which w(p e ) 7^ and w'(p e ) = 0, we have the expansions 

1 w" 

V! = ~^) 2 x 3 + 0(x 5 ), w = w e + 0(x 2 ), 

a = a e + 0{x 2 ) , b = b e + 0{x 2 ) , (19) 

where w" e = e 2be V, w (w e ) / \AR 2 ). The parameters w e and R e must fulfill the constraint equa- 
tion 

3 Pulsation Equations 

We now linearize the field equations on a static background; that is, we write 

w(t,p) = w{p) + Sw(t,p), (20) 
p(t,p) = p(p) + Sn(t,p), (21) 
9ab(t,p) = g ab (p) + Sg ab (t,p), (22) 
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and require that w(p), /x(p) and g a b(p) solve the static equations (|T^)-([T6|). Here and in the 
following, first order quantities are highlighted by boldface letters. 

For solutions without equator, the function R is an admissible coordinate and, moreover, 
one may choose the gauge <5/z = to analyze stability properties. The system of linearized 
field equations then reduces to one pulsation equation (the linearized YM equation) and a set 
of constraint equations (the linearized Einstein equations) . By virtue of these constraints one 
can eliminate the gravitational perturbations from the YM equation, which yields a pulsation 
equation for 5w alone. This technique was used to investigate the stability properties of the 
Bartnik-McKinnon solitons and the SU(2) EYM black holes with hair [0, 0. A detailed 
account of the procedure for a general class of matter models was recently presented in M . 

The Schwarzschild coordinate R is not suited to describe spatially compact solutions. In 
addition, the gauge 5/j, = is not regular at the equator. The aim of the following analysis 
is to generalize the procedure mentioned above, without adopting the gauge dfj, = 0. We 
do so by deriving a pulsation equation for a gauge invariant quantity say, which basically 
reduces to 5w in the Schwarzschild gauge, Sfi = 0. We then argue that there exists an 
everywhere regular gauge for which the metric and matter perturbations can be constructed 
explicitly from £. 

We start by linearizing the relevant field equations in a static background. In all what 
follows, we use 

USw = e- a+b [e a - b 5w')' - e- 2a dib, 
Uw = e~ a+b (e a - b w')' , 
and similarly for all other first order and static background quantities, respectively. 



3.1 The Yang-Mills Equation 

The dynamical YM equation reads 

Uw = 



4 R? 



V» 



(23) 



For the diagonal metric @ we obtain the following expression for the variation of the 
d'Alembertian on a static background: 



5Uw = adw - 2Uw ■ Sb + (dw \ ddc) , 



where 

Variation of the r.h.s. of eq. 



5c = 5a — Sb . 



(24) 



(25) 



gives V )WW I (4R 2 )5w — 2Hw ■ 5/j,, where we have used the 



static equation ([16]) in the second term. The linearized YM equation now assumes the form 



□ " t^I 5w = 2Uw ■ (5b - 5fi) - (dw I d5c) . 

For later use we also write this in terms of the 4-dimensional d'Alembertian, 

a = a + 2(dfx\d ■) , 



(26) 
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and the perturbation dz, 



Sz = — Sw 

R 



A short computation yields ^Hdw = L\Sz — (Rnjs)dz, and thus 



(27) 



dz 



V., 



2R 3 



(5fj, - db) 



R 



(dw | dSc) 



(28) 



3.2 The Constraint Equation 

All first derivatives of the fields with respect to t enter G 00 and T 00 only quadratically. 
The corresponding linearized Einstein equation assumes therefore the form of a constraint 
equation, that is, it contains no time derivatives of the linearized quantities. On the other 
hand, Gqi and T i already are linear expressions in terms of time derivatives, implying that 
the corresponding linearized Einstein equation is a total derivative of a quantity F, say, with 
respect to time. A lengthy computation then shows that the linearized (OO)-equation is the 
total spatial derivative of F . 

This connection between the linearized (00)- and (Ol)-equations is also established di- 
rectly, as a consequence of the Bianchi identity. In order to derive this result, we consider 
the tensor 



E. 



G 



8ttT, 



(29) 

and linearize the contracted Bianchi identity for E^ u . Since the background equations imply 
that only variations of E^ v contribute, we have 







1 



-gSEi 



Hi 



(30) 



Now using the facts that r* ti , r* • and Y % t - vanish in a static spacetime and that r*. 
T l tt = —-gWgttj, we find for the second term in the above expression 



1 n tt n 
2# 9tt,i, 



g tt SE\ - g*8E\ 



Here we have again used the fact that the tensor E\ vanishes identically (off shell) in a static 
spacetime, implying that g u 5E{ = g^SEP^ = 5(g tfl E^) = 5E jt and g ij 5E\ = 5E tj .) In the 
spherically symmetric case we therefore obtain the identity 







-gSE t 



(31) 



Since, as we shall see in a moment, (^—g5E t p ) is already the time-derivative of a quantity 

F, say, 

F = V=9~5E t p 



J t 1 



we can integrate eq. ( |3"T| ) with respect to t. This shows that - up to a integration function 
of p - the linearized (OO)-equation can be obtained from the linearized (Ol)-equation, 



F' = — 



-g 6E* . 
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For a diagonal background metric and a gauge where the variation of the off-diagonal 
part vanishes as well, we find from eq. (0) 



5G 



tp 



R 



(SR)' - V\ t 6R - R' 5T 



i AFP 

pi 



Taking advantage of the fact that the background is static, this yields 



^o 1 



-2 1 

dt 



-(a+b) 



(<5// + (p — a')Sfj, — ji 8b) 



(32) 



The variation of the corresponding component of the stress-energy tensor is immediately 
found from eq. (El), 



R 



(33) 



Integrating the linearized (Ol)-equation ( |3"2] ) with respect to t gives (up to a function of p) 



w' 



Sn' + - a')Sfj, - p'Sb + —6z = . 

R 



(34) 



where, as earlier, 6z = j^Sw. (Since the integration function has the form of an inhomo- 
geneity, it can be set equal to zero by an appropriate choice of the initial conditions. Indeed, 
the linearized (OO)-equation shows that eq. (0) is fixed up to a constant.) The arguments 
presented above now imply that the linearized (OO)-equation is the derivative of this expres- 
sion with respect to p. Hence, eq. (0) comprises the information of both, the linearized 
(00)- and (Ol)-equations. 



3.3 Gravitational Pulsation Equations 

It remains to linearize the (ll)-equation and the (Ai?)-equations. Since we have already 
exhausted the information from the (0a) components, and since (M, g) is spherically symmet- 
ric, it is sufficient to consider the partial traces g ab G a b = g ab G a b and g AB GAB = R~ 2 § AB Gab, 
and to linearize the resulting equations 

Ufi - (dp | dp) - = -A + ^M. (36) 

(Here we have also used d/i = R~ 1 aR + R~ 2 (dR\dR).) With respect to the diagonal 
parametrization (|5D of g we find 

Sap = adu - 2ap- db + (dp \ d[Sc + 28/j,)) 

and 

~6U = adb - 2aa ■ 6b + e- (a+b) (e a - b Sc')' + (da | dSc) . 
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Also noting that 5(dfi\dp,) = — 28b(dfj,\d/j l ) + 2(d/j l \d8fj,) and using the static equation (|15| ) 
in the form Da — (dfi\dfi) = — -^2, we find the linearized equations 

— 2 (dfi | d8^i) — (dfi | dSc) , (37) 

and 

n(6» + 6b) = -2^ + 2(^-^)56+^ 

+ 2d^(d/i5b) + (dfj,-da\dSc) - A8c. (38) 

Here S denotes the co-derivative with respect to the spacetime metric (for instance, d)(d[i8b) 
= D/i- 56+ (dfi\d8b)) , and A is the 4-dimensional d'Alembertian without the time-derivative 
part, 

Adc = -^(^jdc'y. 

It is worth noticing that the operators (d/i | d-) and —d^dfi-) are formally adjoint to each 
other, that is, for arbitrary functions g and h, we have 

h V dll g - gV j dfM h = d ] (g h dfi) , 

where 

£>d M g = (dfi\ dg) , 

V^g = -rfidfig) = -[{Ufj) + V dfJ ] g . (39) 

3.4 Symmetric Form of the Pulsation Equations 

We shall now write the system of linearized equations (p8|) , fl34|), (|37|), (|38|) in the form 

TD<5« = M5v, (40) 

where T and TVf are 4x4 matrices, M is formally self-adjoint, and 8v parametrizes the 
perturbations - suitably arranged in a 4-vector (see below). Using the operators defined 
above, the spatial derivative of the constraint equation (R41) can be written as 



= - 2>k 8b + Vl_ da 5 ii - A 8fi + [ip^Jt 5* , (41) 

where [j^Pd^ 5z = — d)(^dw 8z). 

The linearized equations (pETP, (|3T|) , fl3"B"| ) and (|2"8"| ) can also be obtained from variations 
of the effective action (expanded to second order in the fluctuations) with respect to 8c = 
8a — 8b, 8b, 8/j, and 8z = R^^^Sw. This implies that 

8v = (8c, 8b, 5/x, 8z) . (42) 
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The system fl41|), (|37|) , (|38|) , (|28|) has now the desired form (|40| ), with 

M = P + D , 



(43) 



where P is the symmetric potential matrix 
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— RH-5 = — (dfi\dfi)), and is the formally self-adjoint differential operator 
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The matrix T in front of the d'Alembertian is 



f 
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1 
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1 





V o 








-1 / 



Since the system contains one constraint equation, the rank of T is three. In addition, —T 
has only two positive eigenvalues, reflecting the fact that the effective action for the system 
is indefinite. In fact, the kinetic term in the effective action is 



R 2 e b - a 



R 2 



(dw) 2 + (6b) 2 - (dfi + Sbf 



Hence, the pulsation equations do not form a hyperbolic system in standard representation. 
It is therefore not clear how to make sense of the notion of stability on the basis of the above 
equations. We shall now present a way out of this difficulty, by arguing that the linearized 
system contains actually only one relevant dynamical degree of freedom. The aim is to find a 
gauge invariant generalization of the elimination procedure used in the Schwarzschild gauge 
(<5/x = 0). 



4 Gauge Invariant Formulation 

4.1 Gauge Transformations 

Until now we have not fixed a gauge for the spacetime metric (apart from the vanishing of 
the shift and its variation). In order to isolate the physical degrees of freedom it is, however, 
crucial to get rid of the gauge modes. The metric perturbations and the scalar YM amplitude 
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Sw transform according to (Lxg)^ and Lxw. For a static background with diagonal metric, 
this implies 

8a — > 8a + a f + g , 
8b — ► 8b +b'f + f , 
Sfj, — > Sfi + fi' f , 

Sw — ► Sw +w'f, (44) 

with / = X p and g = X*. Since we require that the shift and its perturbation vanish, / and 
g are subject to 

e 2b f = e 2a g>. (45) 

It is important to note the following: The perturbation da enters the linearized equations 
only via the combination 8a — 8b' = 8c . Using eq. (45), one observes that the gauge 
transformation for this quantity does not involve the gauge freedom g: 



8c' 



8c + 2(a' - b')f + {a" - b")f - e 2b Uf . 



(46) 



Hence, the gauge transformations of the relevant quantities - i.e., the perturbations 8b, 8fi, 
8w and 8c which enter the field equations - involve only one degree of freedom, /. 



4.2 The Schwarzschild Gauge 

As a first application, we consider the gauge 8R = 0, i.e., 8/j, = 0, which we shall also call 
the Schwarzschild gauge. The gauge transformation which achieves 8fi = is regular as 
long as fi' 7^ 0, that is, as long as R is an admissible coordinate for the background solution. 
An advantage of this gauge is that the linearized gravitational equation fl37D becomes also a 
constraint equation. We therefore end up with the YM equation, 



□ 



4R 2 



8w = 2 Ow ■ 8b 



-2b 



W dc 



and the constraint equations ( 34|) and (|3 



yl 8b 

~ 2b y'8c' 



R 2 



w' Sw . 
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--^ aw ■ 8w 
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2a/i-8b. 



(47) 
(48) 



Here we have also used the static background YM equation. Hence, as is well known (see, 
e.g., 0), the gravitational perturbations can be eliminated in the gauge <5/x = 0. The 
linearized YM equation therefore assumes the form of a Schrodinger equation for 8w, 



a 



v., 



AR 2 



8w 



2w' 
y' 2 R 2 



2/i'aw - w'Uy 



• Sw . 



Now using e 2b D/i = fi" + {b' - a' + 2//)// and e 2b aw 
equation for Sw can be brought into the simple form 



w" + {b' — a')w', this pulsation 



□ Sw 



^IWW 

Jr 2 



+ 2e~ (a+6) 



>-v 2 Y 



RR' 



Sw . 



(49) 
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If R can be chosen as a coordinate, then this equation is suited to discuss the perturbations 
of the system. However, if, for instance, the static hypersurfaces are topological 3-spheres, 
then R' has a zero, where the second potential term in eq. (|49|) becomes unbounded (see 
below). 



4.3 The Gauge Invariant Pulsation Equation 



In the previous section we have seen that - in a gauge where dfi vanishes - the system of 
linearized equations can be reduced to a single pulsation equation. We shall now argue that 
this can always be achieved. More precisely, we show that one can find a pulsation equation 
which involves only a particular, gauge invariant combination of the perturbations. It will 
turn out to be convenient to write eq. (fE3) as an equation for fi'dw: 



□ 



(fi'Sw) 



4 R 2 



+ 2e- {a+b) 



RR' 



It is immediately observed from the transformation laws 

£ = n' dw — w' d/i 



,- 1 
i that 



(fx'Sw). (50) 



(51) 



is a gauge invariant combination of the perturbations. Since £ and fi'dw are identical in 
the Schwarzschild gauge (S/j, = 0), and since eq. (|50"D was derived within this gauge, we 
conclude that eq. (|50"D is the pulsation equation for £. 

Using the static equation flT3|), R" / R! = a' + b' — w' 2 /(RR'), the potential on the r.h.s. 
of the pulsation equation ( j50|) can also be cast in the form 



V., 



AR 2 



+ 2D(a + 6-ln| J R , |) - //□— 



Combining the second and third terms in this bracket yields the final result 



AR 2 



+ 2 □ (a + b - n) 



Ofj,' 



•c- 



(52) 



Before we discuss this equation, let us note the following: It is not hard to see that the 
transformation laws (H3) imply that 



fi' 2 db + (// - it'll) dfi - y! dfj,' 



(53) 



is gauge invariant as well. However, by virtue of the constraint equation (p4|) and the 
static background equation (|T3"D, one easily sees that this is equal to R~ 2 w'C- In fact, the 
circumstance that we are able to reconstruct all perturbations from £, and the fact that the 
gauge transformations for the relevant quantities contain only one free function, /, strongly 
suggest that the fluctuation problem involves exactly one relevant gauge invariant quantity. 
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4.4 The Sturm— Liouville Equation 



Let us now consider the pulsation equation fl52| ) for the gauge invariant quantity As 
expected, this is a Sturm-Liouville equation. Indeed, if we introduce the operator 



V 



,a— b 



-d n 



+ q{p) 



with potential 



q(p) = —a 



v. 



□p' 



then eq. becomes 



C + AC = o. 



(54) 

(55) 
(56) 



Here we would like to comment on a subtlety arising in the study of this equation. A 
standard approach is to first consider eq. (|56|) as an ordinary differential equation in some 
L 2 -space. In a second step one then discusses under which conditions a Hilbert space solution 
can be identified with a strict solution in the ordinary sense. In our case, the differential 
equation corresponding to the operator A obviously has singularities at the origin (p = 0) 
and at the equator (p = p e ). (In fact, both points are regular singular points; see below.) It 
is therefore natural to consider A as an operator on the Hilbert space L 2 (7, dr/), where / is 
the open interval (0, p e ) and drj denotes the measure e a ~ b fi'~ 2 dp. Then A is symmetric on 
a dense domain and, as we will show shortly, for both boundary points the limit point case 
(in the sense of Weyl) occurs. Now Weyl's alternative (see, e.g., ||) implies that there is 
precisely one self-adjoint realization of A. Moreover, this realization is essentially self-adjoint 



on 



v = {ue C\l) n L 2 (/, dfi) | Au e L 2 (/, dn) } 



In particular, all functions in T> vanish at the origin and at the equator. On the other 
hand, as we shall demonstrate later, there exist also analytic solutions with C(p e ) 7^ which 
give rise to physically acceptable metric perturbations. By adopting the functional analytic 
approach, we would thus loose half of the interesting modes. 

In order to establish that the equator is a regular singular point, we note that the factors 
in front of and £ diverge not stronger than x~ l and x~ 2 , respectively. More precisely, we 
have 



c" - 



fjr_ 
p' 



+ 0(x) 



p' 



+ 0{l) 



(57) 



where p' = 0(x n ), with n = 1 for odd and n 
indicial polynomial for the above equation is 



3 for even background configurations. The 



r (r — 1) — 2rn + n{n — 1) 



which has the roots 



1 + 2n ± VTT 



1/? 



3, odd case, 
6, 1 even case. 
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Together with the symmetry of the background solution, this implies that the fundamental 
system of eq. fl52|) is analytic near the equator. For r± — 3, this has an expansion of the 
form 



j=0 



j=0 



with cq = ko — 1. For r± = 6, 1 one has 



j=0 



2.7 



2.7 



(5? 



(59) 



with Co = = 1- 

Let us also discuss the origin. To leading order, eq. (|52|) reads 



C" + 



2- + 0(1) 
P 



- + °(p 

p 2 



-1> 



(60) 



With a similar reasoning as above, we now obtain the fundamental solutions 

oo 

C+ = P E c i ^ ' C- 

i=o 



-| oo 
P j=0 



(61) 



It is now easy to verify that both boundary points of the operator A belong to the limit 
point case. For p = p e one has to show that the differential equation Au = has a solution 
for which 

/ 6 \u\ 2 dr] = oo (62) 



for some e > (see, e.g., [|]). (For p = a solution with the corresponding property in the 
vicinity of the origin has to be found.) Using the behavior of the measure dr] = e b ~ a p'~~ 2 dp at 
the singular points and the above expansions for the fundamental systems establishes that 
for both p = and p = p e , A is in the limit point case. 



5 Reconstruction of the Perturbations 

It remains to show that every solution of the pulsation equation for £ gives rise to a regular 
set (8w,8fi,8b,dc) of perturbations. One may try to achieve this in a gauge where the 
relevant equations reduce to a system of algebraic equations for (Sw, dfi,db, 8c) in terms of 
£. The two gauges of this kind are the Schwarzschild gauge, 8/j, = 0, and the gauge 8w = 0. 
However, none of these gauges yields a globally regular description for solutions which posess 
an equator, which demonstrates that this method does not always work. 

A better way is to make use of a suitable gauge for which the perturbations can be 
obtained from a system of ordinary differential equations. This is the case for any gauge 
of the form 8b = constant ■ 8fi, because one can then eliminate the d'Alembertians of 8/j, 
and 8b form the pulsation equations (|3"T| ) and (p8|). In this way one obtains an additional 
ordinary differential equation for 8c. There are two representatives within this class which 
yield globally regular expressions for the perturbations in terms of £. These are 8b = and 
the conformal gauge, 8b = 8fi. 
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5.1 Coverings with two Gauges 



Before we shall work in the conformal gauge, let us briefly consider the two gauges for 
which the perturbations are algebraic expressions in terms of Off the equator, //' does 
not vanish and the Schwarzschild gauge, <5/x = 0, is therefore perfectly regular. (The gauge 
transformation 8/j, — > = 8[i + /// is regular for fjf ^ 0.) Using eqs. fl47|) and fl48|) one 
easily finds for the perturbations in terms of C, 

Sfj, — , = — - , (63) 

Sb=^6 W , ^=-2(^)'^. (64) 

If the background solution has no equator, then the Schwarzschild gauge is everywhere 
regular and the above formulae can be used to compute the perturbations in the entire 
spacetime. If, however, we are considering solutions where R assumes a local maximum, we 
need another gauge in the vicinity of the equator. The odd case, w' e ^ 0, is very simple to 
handle, since one can then choose the gauge 8w = and express 8fi in terms of the gauge 
invariant quantity £• Then using the constraint equation (|34]), we obtain 8b in terms of 
Eventually, taking advantage of the linearized YM equation ( p6|) in the gauge 8w = 0, we 
also obtain 8c' in terms of 

dw = 0, Sfj, = - — , (65) 

w' 

Sb = u = 2e2t J^(£2M. (66) 

Re- a n' w'e~ a // v ' 

It is clear that dfi remains finite since, by assumption, w' does not vanish in a sufficiently 
small neighborhood of the equator. Moreover, we have 8fi = constant + 0(x 2 ), since both 
£ and w' are even. Since the background function Re~ a is even as well, we conclude that 
the numerators in the expressions for 8b and 8c' are of 0(x). Also using // = 0(x) shows 
that both 8b and 8c are well behaved in the vicinity of the equator. 



5.2 A Globally Regular Gauge 

The gauge 8 w = is not suited to compute the perturbations in the even case, w' e = 0, since 
the expressions in eq. (|66|) diverge if // = 0(x 3 ). We now show that there exists a regular 
gauge which gives rise to analytic expressions in the vicinity of the equator for both odd and 
even background configurations. Moreover, this gauge is also regular at the poles, where it 
yields analytic expressions for the perturbations as well. It is therefore possible to cover the 
entire spacetime using just one gauge. 

We consider conformal perturbations of the spatial part of the metric, that is - without 
loss of generality - we choose the gauge 

8b = 8fi. (67) 

(We also note that the gauge transformation between a globally regular gauge and the 
conformal gauge is everywhere regular.) 
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We have already mentioned that - by virtue of any algebraic relation between 8b and 
8 fi - it is possible to eliminate the d'Alembertians from the pulsation equations fl3T|) and 
(j38[), which therefore provides an additional ordinary differential equation for 8c'. For a 
given solution £ of the gauge invariant pulsation equation one can then reconstruct all 
perturbations from £ by solving a system of constraint equations. The relevant equations 
are the pulsation equation (|52| ) for £, the algebraic relations ( |6"TD and (0), the constraint 
([34]) and the new constraint (|7l|) , obtained from eqs. (|37|) and (p8|) : 



□ _ 2e- 2b ^-d p 



AR 2 



2 □ (a + b - ji) - 



[1 dw — w' Sfl = C 



8fl' 

„2a-b 



W' 



a' 8 ii + ^— 8w 



R 



8c' 



W8w + M8fi 



(68) 

(69) 
(70) 

(71) 



where 



W 
M 



6 R- 3 e 2a+b \dw 



-2R- l e 2a+b 



(dw\djj,)} , 
[□(// + 2a) - 3(da\dfi) 



-1] 
R 2 - 



It is worth noticing that the gauge transformation / which achieves 8b 
everywhere regular, in particular at the equator (R! = 0) and at the poles (R 
one easily finds from eqs. 



f = -Re~ b J R^^iSb-Sfi). 



(72) 

<5/x — ^ is 

0). In fact, 



(73) 



The regularity at the equator is manifest, whereas the analyticity at the poles follows from the 
fact that 8b — <5/z = 0(R). (Use the Schwarzschild gauge <5/x = for which 8b = O(R) and 
the fact that the transformation from the Schwarzschild to the conformal gauge is 0(R 2 ).) 

^From a given solution £ of the gauge invariant pulsation equation we now obtain 8 \i 
and 8w from the definition ( |5"9"D of £ and the constraint ([70|): 

<-p e b w' 



8fx 



8n(p ,t)-R>e- b J - — C 



8w = 8w(p ,t) + w' Re b [ 



po R R' 2 

p ga+fc 

po R' 



l 



w 



(74) 
(75) 



where we have also used the background equation (Tl3|). Integrating the additional constraint 
equation (]71f) finally gives 8c': 



8c = 8c'(p ,t) + Re 



b~2a 



[W8w + M8fi] 



(76) 



The time dependent quantities Sfi(po,t), 8w(p Q ,t) and 8c'(p ,t) are not independent. The 
fact that they are subject to the YM equation PS|), together with the residual gauge freedom, 
implies that they are uniquely determined by £ (see the next section). 
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The perturbations obtained from eqs. (74])-(l76D are analytic in the vicinity of the equator 
for both odd and even background configurations: In the even case, where w' = 0(x), R' = 
0(x 3 ) and C — @( x ), the integrands in eq. (|74|) and (|75|) diverge like 0(x~ 4 ) and 0(x~ 2 ), 
respectively Hence, Sfi = 0(x 3 )/0(x~ 4 ) = 0(1) and 6w = O(x) J 0(x~ 2 ) = 0(1). In the 
odd case, where w' = 0(1), R' = O(x) and C = 0(1), one has <5/x = 0(x) J 0{x- 2 ) = 0(1) 
and (5™ = 0(1)/ 0(1) = 0(x). 

The above formulae give also rise to analytic expressions in the vicinity of the poles. 
Using the behavior of the background solution for R — > we have already concluded from 
eq. (H) that C = 0(R). Now using w' = 0(R) shows that <5/x = 0(1) and 6w = O(R). 
Since the terms of 0(R~ 3 ) in M cancel, we have Mdfi = 0(R~ 2 ) and Wdw = 0(R~ 2 ). 
Hence, 8c = 0(R) / 0(^) = 0(1). 

In conclusion, we have shown that the gauge 6b = 8fi gives everywhere regular expres- 
sions for the entire set of perturbations obtained from the gauge invariant quantity £. The 
same also holds for the gauge db = 0, which turns out to be particularly suited for the 
numerical investigation (see Sect. 5.4). 



5.3 The Residual Gauge Freedom 



It remains to discuss the implications of the residual gauge freedom. The transformation 
laws fl44D imply 

(Sb-Sfi) ^(Sb-Sfi) + (b'-fi')f+ /'. 

The conformal gauge Sb — S fi = is therefore invariant under transformations of the form 
f(p,t) = Re~ b <j)(t), where <fi{t) is an arbitrary function which does not depend on p. The 
remaining gauge freedom is thus 



Sw 

8c' 



5v +R'e- b (j){t) , 
Sw +w' Re~ b (p(t) , 

8c + Re b - 2a '(j){t) ~ 



3 a - b (Re- a )' 



dc + Re b - 2a [(f) - k(p) ■ 0] 



0(t) 



(77) 



with 



k(p) 



e a ~ b (Re- a ) 



(n 



Comparing these expressions with the integral formulas (|74])-(|76D shows that contribu- 
tions to dfi of the form R'e~ b a(t) can be compensated by a gauge transformation with 
<fi{t) = —a(t). (The same holds for terms of the form w'Re~ b f3(t) in the expression for Sw.) 
Hence, by virtue of the residual gauge freedom, dfi is uniquely determined by £. Since we 
already know from the integral expressions that 8w and 5c' are regular, we can compute 
these perturbations from the definition ( [5"TD of £ and the linearized YM equation. This 
shows that - once the ambiguity in Sfi has been removed by the residual gauge freedom - 
all perturbations are uniquely determined in terms of £. 
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Figure 1: Bound state energy versus cosmo logical constant for the n = 1 EYM solutions. 



5.4 Numerical Results 



As already mentioned, the gauge Sb = is particularly suited for the numerical analysis. 
Assuming a harmonic time-dependence e lat for the fluctuations, the YM equation (^) and 
the constraint ([34]) become (with 8c = 5a) 



—Q— - V,ww 

dp 4 dp 4R 2 



5w + (Qw'YSfj, + (Qw')Sa' = -Q- l a 2 6w 



(79) 



w 



0*' 



<5/x = 0. 



Here we have also used the background gauge e 2a 
bertians from eqs. 



e = Q(p). Eliminating the d'Alem- 



and (|38f) , we obtain the additional constraint equation 



Q-^(Q^8a!)' ~^Sw- 2e~ 2b ^'Sfi' + ^-(2V-R 2 )d» = 0. 



i? 4 



The main advantage of the gauge 6b = is that the above equations are manifestly regular 
at the critical points R' — 0. They can therefore be used to determine the spectrum for any 
background solution with n nodes and cosmological constant A(n) € [0, A 



rog 



The numerical analysis shows that there exists exactly one boundstate (in the gravita- 
tional sector under consideration) for each solution with n = 1 and A(l) e [0, A reg (l)]. The 
energy of the unstable mode versus the cosmological constant is shown in Fig. 1. For A = 
we obtain the boundstate belonging to the lowest (n = 1) BK solution (in the normalization 
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Q(0) — 1, K — 1). In the other limiting case, that is for A(l) = 3/2, we find a 2 = —2, which 
is the eigenvalue belonging to the unstable mode of the lowest compact solution (see also 
the next section). 

A careful numerical investigation for the few lowest solutions strongly suggests that there 
exist n unstable modes (in the gravitational sector) for any background solution with n zeroes 
of w. 



6 Instability of the n = 1 Compact Solution 

6.1 Gauge Invariant Approach 

The stability analysis of the ground state solution, 

a = b = , i? = sin p , w = cos p , (82) 

can be performed analytically. Using □£ = — £ + Hp = p" = — l/sin 2 p, ^rQp' = 
p"'/p' = 2/sin 2 p and (V, ww )/(AR 2 ) = (3cos 2 p - l)/sin 2 p, one finds from eq. (ID 

-i + < r + ^—c = ^=^<. (83) 

sm p cos p sm p 

The solutions of this equation of the form C(^?p) = C(fc)(p) el<Jfe * are 

C W W=cosV(^^Y, (84) 
v ' \sm p cos p J 

with the spectrum 

a 2 k = k 2 - 3 , where k = 1, 2, . . . . (85) 
In the conformal gauge, 8b = S/j,, we now obtain S/j, from eq. ([74]), 

rp 1 sin(A;p) 
cos 2 p sin p 



c /" p 1 > sm(kp) 

on = cosp / — — C = h cosp-am, (86) 

,/ cos^jO sino 



where a;(t) is a free function. In order to obtain 6w, we can either use eq. (|75|) or the 
definition of C or - most easily - the constraint ([70|), <5iu = sin pdfj,': 

5w = S inp f _ sin 2 . a(t) . (87) 

\ smp J 

Finally, we have to evaluate the integral ( |76]) for dd . This is, however, trivial since the 
background quantities W and M. defined in eq. (|72| ) vanish for the ground state solution 



W=^- 3 (w"-w'p') = 0, M = ~0jr + ^) = 0. 
We therefore conclude from eq. ([76]) that (5c'/ sinp)' = 0, that is, 

Sc = sinp • 7(£) . 
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It is clear that j(t) is not an independent free function, since dd is uniquely determined 
from 8w via the linearized YM equation (|26|). Using the above solution for 8w we find 



□ 



which shows that 



AR 2 



sin p ■ a) 



a — a[ 



d 2 3 cos 2 p — 1 , 
dp 2 



sin 2 p 



sin 2 p 



(d + a) sin 2 p = —w' dd = sin p dd 



7 



d + a . 



(89) 

Eventually, we can use the remaining gauge freedom to get rid of the free function a(t): 
Considering transformations of the form / = Re~ b <p(t) = sin p- <p(t) with <p(t) = —a(t) yields 



Sw 

dd 



Sfi — cos p • a(t) , 
8w + sin 2 p • a(t) , 
Sc — sin p ■ (d + a) 



Sc — sinp • 7(£) . 



(90) 



This implies that - for the ground state - it is consistent to consider conformal perturbations 
of the entire spacetime metric 



Sa = db 



sin( kp) ■ . 
Sn = } HJ e laht 



sin p 

Sw = 8n' sin p e takt 



(91) 



It is also worth noticing that the residual gauge freedom reflects the existence of the conformal 
Killing fields of S 3 . 

The lowest mode, k = 1, corresponds to the negative eigenvalue a 2 ^ = —2 and gives 
rise to exponentially growing perturbations. As is seen from the above solution, the metric 
perturbations depend only on the time coordinate and the YM amplitude is not perturbed 
at all for the unstable mode, 



Sa = db = 6n = e Vlt , dw = 



(92) 



The second lowest mode, i.e. the mode with a 2 2 ^ 



c 



(fc=2) 



cos 2 p( 



sin(2p) > 
sin p cos p 1 



+1, is a pure gauge mode, since 
0. In fact, we find from eqs. (^6[), ( ^7[) and (|88D for k = 2: 



dfi = 2 cos p e* 



Sw 



-2 sin 2 p e u , 8c = . 



A gauge transformation with / = —2 sin pe lt then gives 6/j, = dw = 8c = 0. 



6.2 An Alternative Approach 

The spectrum of perturbations of the ground state is also obtained without using the gauge 
invariant quantity £ as follows: Since the spatial part of the background solution is a space 
of constant curvature (the geometrical 3-sphere), it is natural to introduce the two gauge 
invariant Bardeen potentials used in cosmology to describe scalar perturbations |12[ (see, 
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e.g., for a review). The Bardeen potentials are proportional to 8a and 8b = dfi. We are 
therefore again in the conformal gauge. The fact that the norm of the static Killing field of 
the background solution is constant (a = 0) implies that the remaining gauge freedom can 
be used to set 8c = 0, i.e. 8a = 8b. (This is a peculiarity of the ground state, for which the 
equation for 8c' decouples and implies that 6c is actually a pure gauge; see above.) Hence, 
by virtue of the special form of the background solution, there exists only one independent 
scalar Bardeen potential. It is therefore consistent to consider perturbations of the form 

8a = Sb = dfj, , Sfj,' = 8z, (93) 

where the last equation is the constraint (^) in the gauge 8b = 8fi (with a' = 0, R = —w' = 
sin p.) 

Due to the conformal invariance in four dimensions, the metric perturbations do not 
enter the linearized YM equation (p8|) for 8a = 8b = 8/j,, 



cos 2 p 
sin p 



8z = 0. (94) 



By virtue of the relations (|9~3"D, we obtain the following pulsation equation for 8n from eq. 

m 

[□ + 2]Sfj, = 0. (95) 

Now using (ndfj.)' = nS/j! + 2(ff£) / V = (□ - ^r)8fj! and 8^' = 8z, shows that the 
pulsation equation for 8w is the derivative of the pulsation equation for 8fi: 



([□ + 2] 5rf= (o- 2 g)fc.O. 

Since the background solution is a geometrical 3-sphere, the 4-dimensional d'Alembertian 
becomes □ = sin _2 p<9p(sin 2 p<9 p ). The solutions of eq. ( |95|) are therefore the harmonics of 

S 3 , 

sin((l+l)p) e „ 
sin p 

where the spectrum is shifted by —2: 

a\ = I (£ + 2) - 2 = [I + l) 2 - 3 . (97) 
Since the unstable mode {I = 0) is purely time dependent, 8fi = e^ 2 *, and since the YM 



equation ( 94]) is the derivative of eq. (p5|), the unstable mode is missing in the spectrum of 
the YM equation. Using 8w/ sin p = 8z = 8fi', p! = cos pj sin p and w' = — sin p, we finally 
find 

£ = p'8w — w'Sfi = cos p ■ 8fi' + sin p ■ 8fi 



,2 



8» V 2 /sin((£+l)p)V 



= cos pi — — I = cos p . . . 

ycos^py y snip cos p J 

(98) 

which - with k = i + 1 - is in agreement with the results presented in the previous section. 
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6.3 SUSY Transformation 



The potential appearing in the pulsation equation (f49|) for dw (or in the corresponding 
equation ( |5"2"D for £) can diverge at the equator. (We have argued earlier that the equator 
is a regular singular point.) A pulsation equation with bounded potential can be obtained 
after a supersymmetric transformation. This requires, however, the knowledge of at least 
one solution of the original equation. Unfortunately, we were not able to find such a special 
solution in the general case. Nevertheless, the method works perfectly for the ground state. 

In order to demonstrate this, we consider the pulsation equation fl49|) for dw. Since this 
equation was derived in the Schwarzschild gauge, dfj, = 0, and since 



1 . w' 
— Q = 5w jdfi 

ft' pJ 



(99) 



is clearly gauge invariant, eq. (HDp is also the correct equation for £. For the S 3 background 
solution we have □ = —d 2 /dt 2 + d 2 /dp 2 , and therefore 



_a 2 

dp 2 



(100) 



where the potential diverges at the equator, p e = ir/2, 



V 



3 cos 2 p — 1 
sin 2 p 



+ 



cos 2 p 



(101) 



Let if) be a particular solution of eq. ( |100|) with eigenvalue a 2 , say. Then the differential 
operators B^ and B^ are defined as 



* dp V ' 



(102) 



Now using V = <ri + ip f/ /ip and BtB^ = d 2 /dp 2 — ip"/ip, eq. (|100|) assumes the form 



B% B- 4 , C 



cr^ - a 



Multiplying with B^ gives 



where r)^, is the supersymmetric partner of C, with respect to if), 

■7 t> $ > 



Ip 



'Hip — XJ ijj 

Since B^B+ = d 2 /dp 2 + (ip'/if))' - (*p'/ip) 2 , eq. (|TU§ finally becomes 



(103) 



(104) 



(105) 



dp 2 



- V + 2 



'if 



lip 



v 2 ^ 



(106) 
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Hence, the potential in eq. (|100|) picks up the additional contribution — 2(?///?/>) / which can 
compensate the divergent term in V. (It is obvious from the above derivation that the 
eigenvalue a\ does not lie in the spectrum of the equation for the supersymmetric partner 
with respect to ip.) 

As an example, we consider for ip the lowest mode, ip — Cm = C/V = sin 2 p/ cos p. The 
effective potential for rj^ then becomes finite at the equator, 

\ if) J sin p 

Since we have performed a supersymmetric transformation with respect to C,{i)i the unstable 
mode is missing in eq. (|106|) with potential ( 107 ). In order to obtain this mode from an 
equation with finite potential, we also consider the transformation with respect to the third 
mode, £(3) (recall that 2 ) = 0): 

I . ( sin(3p) V (c 2 -!)(! + 4c 2 ) 

v = — C(3) = sm P cos p = > 

p' w y sin p cos p J c 

where c = cos p. A straightforward calculation now yields the supersymmetric potential 

(f\ 48c 6 - 24c 4 + 139c 2 - 13 
V ~ 2 U) = (l-c 2 )(l + 4c 2 ) 2 ' (1 ° 8) 



which is bounded in the vicinity of the equator. The pulsation equation ( |106|) with potential 
08|) has the unstable mode 



8 sin 3 p 
r 1 + 4 cos^p 

with negative eigenvalue a 2 = —2. (It is now an easy task to verify that 77^ is indeed the 
supersymmetric partner of C(i)/A i ' = sin 2 p/ cosp with respect to ip, 



's 2 \ (f(i+4gy s 2 



f(l+4c 2 ) 



where s = sin p.) Hence, we have obtained the unstable mode from a perfectly regular 
pulsation equation in standard form. 



7 The Odd-Parity Sector 

In order to study purely magnetic solutions it is sufficient to work with a gauge potential of 
the form (f|). However, if one is interested in the stability properties of such configurations, 
one must, in general, consider perturbations of the full spherically symmetric gauge potential 

A = At p + zu {f^d'd + Ttp sm.'ddip) + (w — 1) {f Lp d'd — sm.'ddip) , (HO) 
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where A = a^dt + a\dp. It is a crucial observation that 8 A and 8vj decouple from the 
remaining perturbations, that is, the equations for 8 A and 8~cj contain no metric pertur- 
bations and no 8w's - and vice versa. This is a consequence of the fact that the parity 
transformation 

■& -> 7T - ft, if -> 7T + <fi (111) 

acts as a symmetry operation on the background solutions. Until now, we have restricted 
ourselves to the fluctuations with even parity, that is, we have considered perturbations for 
which 8 A = and 8~cj = 0. 

In this section we shall investigate the orthogonal set of fluctuations that is, we discuss 
the pulsation equations for 8A and 8z& on a static, purely magnetic background. We do 
so by taking advantage of a method developed in fLlf . The following analysis applies to all 



solutions of the EYM equations with cosmological constant studied in M. The reason for 
this lies in the fact that it is actually sufficient to consider perturbations which vanish at the 
horizon (see fl5j for details.) 



7.1 Pulsation Equations 



The YM equations for the general spherically symmetric gauge potential ( |110| ) are most easily 
obtained from the YM action, S = f M tr(F A *F). Using the general form g = g + R 2 g of 
the metric and the fact that w, m and R are functions on M, whereas A is a 1-form on M, 
one finds 



S 



where V(w, w) 



M 



2 



dA I 2 + I dm - Aw I 2 + I dw + Am 



2 , V{w,m) 



2R 2 



(m 2 + w 2 



I) 2 and I dA 



2 fj 



dt 



—e 



-2a 



dr 



- dA A *dA, f] = y 7 ^ dt A dr. Now using 
26 , and A = a^dt + a\dp , 



we obtain the YM equations upon performing variations of S with respect to a , ai, m and 
w. The linearization of these equations on a static background with = ai = m = finally 
gives 

-\ i 



R 2 e- {a+b \8d 1 -8a ') 
R 2 e~ {a+b \8d 1 -8a '\ 



-2e 
-2e 



b— a 



a—b 



w 8a — w8tzj 



w 2 8a 1 — w8zu' + w'8zu 



e b - a (8^J - w8a ) 



+ 



e a - h w'8ax + e a+b 



e a - b (8zu' - w8ax) 

2 -l 



w 



R 2 



8zu . 



(112) 
(113) 
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In addition, we have the linearized YM equation (p6[) for 8w, which remains unchanged. As 
we already pointed out, the above three equations for the fluctuations with odd parity do 
neither involve the metric perturbations nor the even parity quantity 8w. 

We now adopt the temporal gauge and consider the following separation ansatz for the 
perturbations: 

8a o = 0, S ai = ^- X (j>)e M , 8^ = i{p)e M . (115) 
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Also introducing the background quantity 7 and the new coordinate x, 



2 = ±_ 2a 

7 i? 2 



dx = e dp, 

yields the following set of ordinary differential equations for x an d 



(116) 







w 2 7 2 X + w, x £ - w£, x = o 2 X- 



(117) 
(118) 

(w 2 7 2 X + w,x £ ~ vj£, x ), x = a 2 w£. (119) 

In the last equation we have also used the background YM equation for w in the form 
2w, xx = 7 2 u>(u> 2 — 1). Equation ( |117| ) is an obvious consequence of the first two equations. 
In order to gain a second order equation from the first order equations (117) and (|118| ), we 
introduce the function cp = xl w i for which we find the Schrodinger equation 



<P, 



1 2/ 2 

-7 (W 

2 ' K 



2 l« 

\ w . 



(p = a 2 <p. 



(120) 



(For cr 2 7^ this equation is equivalent to the system (|117|) -( PT5 ). For cr 2 = we impose eq. 
( |1 2 0|) as a gauge fixing condition; see A particular solution - which is related to the 

residual gauge freedom for a 2 = - is (see |TJ) 



00 



n. 



7 2 w 



:i2v 



where the function Q(x) is subject to the differential equation 



,1 J 



w 2 Vl. 



:i22) 



7.2 The Number of Odd-Parity Instabilities 



The potential in the Schrodinger equation ( |120| ) is non-negative and becomes unbounded 
at the zeroes of the background amplitude w. This suggests that the number of nodes of w 
determines the number of eigenfunctions with negative a 2 . 



In fact, in ITT 



we have shown 

that this is the case for the BK solitons and the corresponding black hole solutions. We 
shall now take advantage of the technique introduced in [TT| to show that the number of 
nodes of w equals the number of unstable modes for all solutions with cosmological constant 
described in [0. 

The key idea is to apply a supersymmetric transformation of eq. (|120|) with respect to the 
zero-energy solution <p given in eq. (|121 ). The supersymmetric partner ip of (p then fulfills 



a similar Schrodinger equatian as (p, where the additional term —2(p , x / <p is added to the 



potential in eq. (|121|) (see also sec. 5.3, eqs. fllOOp , (|106|) ). In |TI| we have demonstrated 



that the pulsation equation for ip can then be written in the form 



"0 )xx H 



-j 2 {3w 2 -l) + 2(w 2 Z), x 



ip = a 2 ip, 



(123) 
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where Z is obtained from 



Z(x) 



X 
ft, 



n + 



C+ l X (w~f 2 /€t, x ) 2 dx 



(124) 



Here C is an arbitrary constant and Q denotes a particular solution of the differential equa- 
tion fll22| ). For <t = 0, eq. (|123|) admits the solution 



i/> (x) 



w exp 



w Z dx 



j-t) 



(125) 



If Z is everywhere regular, then ip is also smooth and regular. This implies that i/> has 
n zeroes, which will enable us to count the number of bound states. Hence, we have to 
investigate the behavior of Q, from which we obtain Z and i/'o by virtue of eqs. ([124]) and 
( |125|) , respectively. 

In order to analyze the above equations, it is convenient to choose the background gauge 



3 2a(p) 



-2b(p) 



Q{P) 



(126) 



which we have also used in [JJ. Also returning to the radial coordinate p (dx = dp/Q, 
7 2 = 2Q/R 2 ), the differential equation ( |122j ) for fl becomes 



2w 2 Q. 



(127) 



Our first aim is to show that - for all static background solutions with horizon discussed 
in |L| - there exists a smooth solution H(p) of eq. ( |127|) which behaves as 



n(p) 



k 

7 



+ 0(1), and Q(p) = 0{p- Ph ) 
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in the vicinity of the origin (p = 0) and the horizon (p = ph), respectively. Moreover, Q 2 is 
monotonically decreasing between p = and p = ph- 

In order to discuss the behavior near the origin, we recall that all solutions obtained in 
behave as (see 0, eq. (49)) 



w 



1 + 0{p 2 ), R = p + 0{p 5 ), Q = 1 + 0{p 2 



(129) 



for p — > 0. Hence, the origin is a regular singular point of eq. ( |127| ) with roots r± = 1, —2. 
In the vicinity of the horizon we have (see |], eq. (51)) 



w 



w h + O(Ap), R = R h + 0(Ap), Q = \Q' h \ Ap + 0(Ap 2 



(130) 



where Ap = ph — p and Q' h < 0. Using this in eq. (|127| ) shows that there exists a local 
solution with Cl(p) = O(Ap). Finally, in order to establish the monotonicity property, we 
multiply eq. ( |117| ) with Q/Q and integrate from p < p h to p. Using fl(p h ) = 0, this yields 
the desired result: 



(n 2 )\ P ) 



R 2 (p) 



R 2 n' 2 + 2^n 2 \dp < o 



(131) 
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Using the expansions for Q and the definition dx = dp/Q of the coordinate x, we find 
x = p + 0(p 3 ) as p — > and x = — \Qh\~ 1 ln(Ap) + 0(1) as p — ► p^. Hence, eq. (|128|) implies 

= 4 + 0(1), and n(x) = 0(exp(-\Q' h \x)) (132) 
x^ 

in the vicinity of the origin (x = 0) and the horizon (x = oo), respectively. Using this in the 



expression ( 124 ) for the function Z(x) shows that 



Z(x) = - + 0(x), and Z(x) = -wl- + 0(^) (133) 

JU JU JU 



as x — > and x — > oo, respectively. Finally, by virtue of eq. ( ]125p , we find for the solution 
ij} Q (x) of eq. O 



if) {x) = c x + 0(x 2 ), and VoW = Coo - + 0{\) (134) 



X ' X 2 ' 



as x — > and x — >• cxd, respectively. This shows that tp G -^ 2 (0, oo). Since ^ is the zero- 
energy solution of eq. ( |123| ), and since ipo has the same number of nodes as the background 
solution w on the interval (0, oo), we conclude that eq. ( |123| ) has n bound states, ipk(x) G 
L 2 (0, oo) with al < 0, k = 1 . . . n. 

For a ^ 0, the supersymmetric relation between (f> and if) (with respect to O ) guarantees 
that eq. ( |120| ) has also n bound states. Finally since the Schrodinger equation for (f) is 
equivalent to the system (|117|) -( |TT9D , the same holds for the solutions x an d ^- Explicitly, 



one finds (see fllTf ) 

Xk = ^ ( Uj3z ' l f ) k + w'^k ~ w ^'k) > £ fc = *l>k ~ wZ X k - (135) 

This proves the existence of exactly n unstable modes in the spherically-symmetric, odd- 
parity sector for all EYM solutions with < A < A*(n) and A*(n) < A < A reg (n) obtained 
in |]T|. (The special case A = A^(n), for which the horizon occurs at the maximum p e of 
R, is not covered by the above reasoning, since eq. ( |127|) develops an irregular singular 
point at p = p e = ph. Although we have not carried out the analysis for this case, we 
expect on continuity grounds that the result remains the same.) In the limit A — > A reg (n), 
corresponding to the compact solutions, the horizon shrinks to a point (the "south pole"), 
where the boundary conditions are identical with the boundary conditions at the origin. The 
above analysis can easily be repeated for this situation, and reveals again n unstable modes 
for the compact solutions with n zeroes. 



8 Concluding Remarks 

In this paper we have investigated the stability of the static, purely magnetic EYM solu- 
tions with regular center and positive cosmological constant. Whilst the numerical analysis 
is straightforward for the asymptotically deSitter solutions, new difficulties arise for the 
compact and the bag of gold type configurations. 
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In order to solve these problems, we have generalized an elimination procedure for the 
gravitational fluctuations - earlier derived for the Schwarzschild gauge (8R = 0) - in a 
gauge invariant manner. In this way, we end up with a standard Sturm-Liouville equation 
for a gauge invariant quantity £• We have finally argued that there exists a globally regular 
gauge, with respect to which the metric and matter fluctuations can be reconstructed from 

c 

The numerical investigation for the first few branches (n = 1,2,3) shows that there 
are n unstable modes for every solution with n zeroes of w (in the even parity sector). 
This turns out to be true for the entire interval [0, A reg (n)], that is, for the BK solutions, 
their asymptotically deSitter-like analogues, the bag of gold solutions and the compact 
configurations. 

Adopting a method originally developed for the BK soliton and black hole solutions, 
we have finally also established analytically the existence of n unstable modes in the odd 
parity sector. The procedure, which is based on a supersymmetric transformation of the 
Schrodinger equation for the sphaleron-like modes, was explicitly carried out for the bag of 
gold solutions. As expected, the compact solutions also exhibit n unstable modes with odd 
parity. Hence, we conclude that every solution - characterized by the number of zeroes n 
and the cosmological constant A(n) G [0, A reg (n)] - has exactly 2n unstable modes. 
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